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NH  {X  32“,  v=  1-3)  formation  and  vibrational  relaxation  in  electron-irradiated 
Ar/N2/H2  mixtures 

James  A.  Dodd,*’  Steven  J.  Upson,  Dorothy  J.  Flanagan, b)  and  William  A.  M.  Blumberg 

Geophysics  Laboratory  (AFSCJ/OPl,  Hanscom  Air  Force  Base.  Massachusetts  01 73 1 

James  C.  Person  and  Byron  David  Green 

Physical  Sciences  Inc.,  20  New  England  Business  Center.  Andover,  Massachusetts  01810 
(Received  24  October  1990;  accepted  14  December  1990) 

Measurements  of  the  dynamics  of  NH(^32  ~,v  =  1-3 ).  created  in  electron-irradiated  N2/H2 
and  Ar/N2/H2  mixtures,  have  been  performed.  Time-resolved  Fourier  spectroscopy  was  used 
to  observe  NH(t>— u-l )  vibrational  fundamental  band  emission.  Time-dependent  populations 
were  then  determined  by  spectral  fitting.  Subsequent  kinetic  fitting  of  these  populations  using  a 
single-quantum  relaxation  model  and  a  power-law  dependence  of  kt,  on  v  yielded  the  following 
NH(u  =  1-3)  relaxation  rate  constants  (units  of  10"  14  cm3  s" '):  ,  (N2)  =  1.2  +  0.5, 

*„.i(N;)  =3.8±  l.5.Jfc„_,(N,)  =7.5±2.5;Ar„_,(Ar)  =0.2  ±0.1, 

2  (Ar)  =  0.5  ±  0.2,  j  ( Ar)  =  0.S  ±  0.3;  ,  (H2)<50.  k„,  2  (H,)<100, 

j  (H-)<  150.  In  addition,  the  N-/H-  data  provided  a  measurement  of  the  nascent  excited 
vibrational  state  distribution  resulting  from  the  reaction  N(:Z?)  -f-  H,  —  NH(Xu)  +  H.  The 
ratio  NH(  1  ):NH(2):NH(3)  was  found  to  be  1.0:0.97:0.81  (  +  0.28  in  each  value). 

Comparison  of  the  observed  nascent  distribution  with  that  of  a  statistical  model  suggests  that 
the  ratio  NH(0):NH(  1 )  =  0.47.  Using  this  derived  distribution,  we  find  the  average  product 
level  ( v )  =  1.6.  and  the  fraction  of  the  available  product  energy  in  vibration  (f„)  —  0.44.  The 
present  evidence  confirms  that  a  single  reaction  mechanism  dominates  NH  formation,  and 
suggests  that  the  reaction  proceeds  by  direct  H  atom  abstraction  rather  than  the  formation  of  a 
long-lived  H-N-H  intermediate. 


I.  INTRODUCTION 

The  imidogen  radical  NH(X'S~^is  a  species  of  funda¬ 
mental  scientific  interest.  NH  is  closely  analogous  to  the 
isoelectronic  triplet  species  CH2  and  ground  state  0(i jP), 
which  have  received  a  great  deal  of  attention.  Also,  because 
NH  has  a  relatively  small  number  of  electrons,  its  electronic 
structure  and  properties  are  amenable  to  calculation.  How¬ 
ever,  aspects  of  the  dynamics  of  NH  remain  poorly  under¬ 
stood.  HPof  Instance,  while  studies  of  excited  state  [e.g., 
NH(a'A)  1  quenchihg^by  buffer  gases  have  been  pub¬ 
lished,1,2  the  rotational  ancf'vihiational  dynamics  of  ground 
state  NH  ( X.v )  have  not  been  reported.  NH  is  known  to  play 
an  important  role  in  the  combjistion  of  nitrogenous  materi¬ 
als,  such  as  hydrazine  fuels."  NH  has  also  been  observed  in 


has  also  been  observed  in 


astrophy$ical  objects,  including  comets  and  certain  types  of 
starsf^Thus,  knowledge  of  NH  reaction  dynamics  is  needed 
to  characterize  NH  in  numerous  chemical  and  physical  pro¬ 
cesses^ 

Studies  have  been  performed  at  the  Geophysics  Labora¬ 
tory  to  understand  the  dynamics  of  formation  and  relaxation 
of  NH(Xo  =  1-3)  created  by  the  interaction  of  a  35  kV 
electron  beam  with  gas  mixtures  containing  N2,  H2,  and  Ar. 
NH(u)  fundamental  band  infrared  emission  was  first  ob¬ 
served  in  electron-irradiated  N2/H2  mixtures  by  Green  and 
Caledonia.’  By  selectively  quenching  pathways  involving 
Jqqic  precursors,  they  determined  that  the  observed  NH(u) 

2  */..»  fu-  rf  .-V7 1  Q )  <  'I 

‘'Currently  at  Stewart  Radiance  Laboratory,  139  Great  Rond.  Bedford. 
MA  01730. 

h|  Current  address:  32  Duggan  Road,  Acton.  MA  01720. 


was  produced  by  the  reaction 

N(:Z?)  +  H2  — NH(u<;3)  +  H.  (1) 

An  N(:Z>)  precursor  is  favored  over  ground  state  N(4S), 
since  N(4S)  +  H2  — NH  +  H  is  endothermic,  and  also  over 
the  higher  energy  N  ( 2F)  state,  for  which  no  adiabatic  path¬ 
way  exists  connecting  reactants  to  products.'’  In  addition, 
the  exothermicity  of  reaction  (1)  permits  formation  of 
NH(u<3),  consistent  with  observation.  Piper  et  al1  have 
measured  the  rate  constant  for  quenching  of  N(2Z>)  by  H2  to 
be  2.3  +  0.5  (  —  12)  em’/s  (read  as  2.3X10-12).  While 
N(3£>)  is  probably  quenched  by  H2  primarily  through  reac¬ 
tion  ( 1  ),7  the  branching  ratio  between  reaction  ( 1 )  and  non¬ 
reactive  quenching  has  not  been  measured. 

We  have  performed  two  series  of  experiments  concern¬ 
ing  NH(u)  formation  and  deactivation,  each  consisting  of 
about  175  runs  under  variable  conditions  of  pressure  and 
electron  beam  parameters.  The  first  series  involves  mixtures 
of  N;  and  H2  only,  while  the  second  involves  mixtures  of  Ar, 
Nj,  and  H:.  Rate  constants  for  the  relaxation  of  NH(y  =  1- 
3)  by  collision  with  Ar,  N2,  and  H2  have  been  determined 
( for  Ar  and  N2)  and  bounded  (for  H2).  In  addition  to  relax¬ 
ation  rate  constants,  the  N2/H2  mixture  data  provide  infor¬ 
mation  on  the  relative  formation  rates  for  NH(u=  1-3). 
Given  the  mechanism  for  formation  of  NH(u)  in  reaction 
( 1 ),  the  observed  relative  formation  rates  can  be  equated  to 
the  nascent  vibrational  distribution.  To  our  knowledge,  this 
experiment  provides  the  first  measurement  for  NH(u)  relax¬ 
ation  rate  constants  and  for  the  nascent  vibrational  distribu¬ 
tion  resulting  from  reaction  ( 1 ). 
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II.  EXPERIMENTAL 

The  experiments  were  performed  using  the  room-tem¬ 
perature  LABCEDE  apparatus,  which  has  been  described  in 
previous  publications.*  ''  A  pulsed  35  kV  electron  beam  ex¬ 
cites  a  mixture  of  gases  in  a  5.4  liter,  stainless  steel  vacuum 
chamber  (“target  chamber"),  connected  to  a  mechanical 
pump.  Computer-produced  pulses  were  transmitted  over  a 
fiber  optic  link  to  a  high-speed  transistor  controlling  the 
electron  beam  grid  voltage,  enabling  beam  switching  times 
of  about  5  fis.  For  the  N2/H2  study,  the  electron  beam  was 
pulsed  on  for  either  300  or  600  /zs,  with  7.5-30  ms  between 
pulses.  The  instantaneous  current  observed  at  the  target 
chamber  was  2.3-3. 3  mA.  For  the  Ar/N2/H2  study,  smaller 
doses  were  used:  pulse  widths  of  300-600  /zs.  repeated  at  20 
ms  intervals,  with  currents  of  0.6-1. 5  mA. 

Argon,  nitrogen,  and  hydrogen  (all  99.999%)  were 
used  without  further  purification.  Calibrated  flow  control¬ 
lers  and  aO-lOOTorr  pressure  transducer  were  used  to  regu¬ 
late  gas  flows  and  partial  pressures.  Flow’  rates  and  pressures 
were  as  follows:  For  the  N2/H:  study.  N2,  3810  seem.  4-25 
Torr;  H2,  20-100  seem,  0.05-0.40  Torr.  For  the  Ar/N2/H2 
study,  Ar,  0-5000  seem,  0-1 8  Torr;  N2. 1 000-5000  seem.  4- 
21  Torr;  H2,  40-500  seem.  0.08-1.00  Torr.  Fast  flow  rates 
and  short  residence  times  were  used  to  minimize  buildup  of 
beam-created  quenchers  in  the  target  chamber.  For  the 
Ar/N2/H:  study  the  gases  were  passed  through  a  coiled  cop¬ 
per  tube  cooled  to  —  78  °C  to  eliminate  possible  H-O  conta¬ 
mination. 

For  N;/H:  mixtures,  observed  fluorescence  both  during 
and  following  beam  excitation  is  dominated  by 
NH(u-i)  —  1 )  and  N,(2?  —  A)  first-positive  band  emission. 
The  N;  first-positive  band  emission  does  not  obscure  NH  ( v ) 
emission.  On  the  other  hand,  electron  irradiation  of 
Ar/N;/H;  mixtures  results  in  intense  Ar*  emission  during 
the  beam  pulse,  which  does  obscure  NH  ( v)  fluorescence:  20 
/zs  after  termination  of  the  beam  pulse,  the  Ar*  fluorescence 
is  fully  quenched,  and  NH(i')  emission  can  be  observed 
without  interference. 

NH(u)  emission  was  observed  at  a  90' angle  to  the  elec¬ 
tron  beam  using  a  slow-scan  Michelson  interferometer  oper¬ 
ated  with  an  InSb  detector.  The  NH(u  — 1>  —  1 )  vibrational 
fundamental  band  emission  region  of  2500-3400  cm"'  is 
near  the  peak  of  the  spectral  response  of  the  detector.  Heli¬ 
um-neon  laser  fringes  were  used  to  monitor  displacement  of 
the  moving  mirror  in  the  interferometer.  Spectra  were  ob¬ 
tained  with  15-20  cm  ~  1  FWHM  resolution,  corresponding 
to  1024  He-Ne  reference  laser  fringes.  Fluorescence  decay 
information  from  50  electron-beam  duty  cycles  was  aver¬ 
aged  during  each  fringe  in  order  to  improve  the  S/N."  For 
the  N:/H;  runs,  the  sample  spacing  was  limited  to  a  mini¬ 
mum  of  25  /zs  by  the  PDP  1 1/40  data-taking  computer;  for 
the  Ar/N2/H;  runs,  a  Masscomp  5300  computer  enabled 
sample  spacings  as  small  as  3  /z s  and  better  characterization 
of  observed  infrared  fluorescence  decay.  The  detector  output 
was  successively  passed  through  a  matched  preamplifier 
(factor  of  10  gain)  and  a  PAR  113  preamplifier  with  1  Hz 
low-frequency  and  10  kHz  high-frequency  rollofT  settings 
(factor  of  200-1000  gain),  and  then  digitized.  After  each 
run,  the  averaged,  time-ordered  interferograms  were  trans¬ 


ferred  to  Apollo  DN3000  and  DN10000  computers  for 
further  analysis. 

Blackbody  radiation  produces  significant  modulation  in 
the  observed  interferograms,  resulting  in  a  spectral  baseline 
which  rises  with  decreasing  frequency  at  wavenum¬ 
bers  <  2500  cm  "  Thus,  to  permit  the  observation  of  possi¬ 
ble  features  in  this  region,  an  interferogram  taken  long  after 
termination  of  the  e  beam  (i.e.,  after  all  molecular  emission 
has  been  quenched)  was  subtracted  from  NH(u)  interfero¬ 
grams.  Fast  Fourier  transformation10  of  the  resultant  inter¬ 
ferograms  generated  real,  time-ordered  fluorescence  spec¬ 
tra.  Each  double-sided  interferogram  was  zero-filled  to  8192 
points  to  give  an  adequate  spectral  point  density.  The  wave- 
number  dependence  of  the  phase  of  the  transform  was  cor¬ 
rected"  by  subtracting  a  reference  phase.  The  reference 
phase  was  determined  by  transforming  a  triangularly  apo- 
dized  center  portion  of  the  interferogram.  corresponding  to 
60  He-Ne  laser  fringes.  Finally,  spectral  fitting  was  per¬ 
formed.  described  in  detail  in  the  next  section. 

Calibration  of  the  spectral  response  of  the  detection  sys¬ 
tem  was  accomplished  using  an  Infrared  Industries  model 
408  source  to  produce  1000  °C  blackbody  radiation.  This 
radiation  was  expanded  with  a  CaF:  telescope  to  fill  the 
ficld-of-view  of  the  detector.  The  response  was  checked  near 
the  end  of  the  study  by  reflecting  the  blackbody  radiation 
with  a  Labsphere  model  IRT-94-020  factory-calibrated 
gold-coated  diffuse  reflector  to  fill  the  field-of-view.  The  two 
methods  gave  consistent  calibration  curves. 

III.  SPECTRAL  FITTING  OF  NH(v) 
CHEMILUMINESCENCE 

A  typical  NH(u)  emission  spectrum  corrected  for  de¬ 
tector  spectral  response  is  shown  in  Fig.  1.  Emission  from 
upper-state  population  in  v  =  1-3  is  evident.  No  fluores¬ 
cence  is  observed  in  the  overtone  region,  which  also  lies  well 
within  the  spectra!  range  of  the  detector.  All  spectra  were  fit 


Wavenumber  (cm-1) 


FIG.  1  Typical  observed  NH(r)  vibrational  fundamental  band  emission 
spectrum,  corrected  for  spectral  response.  Data  points  have  been  connected 
by  straight  tines.  Emission  from  e  =  1-3  is  observed.  A  spectral  fit  is  shown 
offset  beneath  the  spectrum,  zero  tines  are  included  for  both  plots.  Pressures 
used:  21.0  Toir  N-.  0.08  Torr  H..  Electron  beam  parameters:  35.0  kV,  1.5 
mA  at  the  target  chamber.  0.6  ms  duration.  19.5  ms  between  pulses  Spec¬ 
trum  was  taken  30  /rs  after  beam  termination. 
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in  the  wave  number  range  2500-3400  cm " 

The  spectral  analysis  of  NH(u)  is  complicated  by  the 
lack  of  validated,  accurate  Einstein  coefficients  for  infrared 
fundamental  band  emission.  No  calculated  set  is  available  in 
the  literature.  Thus,  we  have  adopted  the  following  ap¬ 
proach:  For  the  relative  vibrational  bandstrengths,  A„  coeffi¬ 
cients  were  obtained  by  computing12  matrix  elements 
|<e|^(r)|u  —  1)|2.  The  dipole  moment  function  /x(r)  of 
Goldfield  and  Kirby413  was  used  for  the  calculation.  Rota¬ 
tionless  wave  functions  (u|  were  computed  in  two  steps. 
First,  an  NH(Af)  potential  energy  function  was  determined 
through  a  Rydberg-Klein-Rees  inversion14  of  a  set  of  em¬ 
pirical  NH(A")  rovibrational  term  values.15  The  wave  func¬ 
tions  (u|  were  then  calculated  from  the  NH(.T)  potential 
function  using  a  Numerov-Cooley  algorithm  to  solve  the 
radial  Schrodinger  equation.1'’  This  procedure  has  been  de¬ 
scribed  in  more  detail  elsewhere  for  the  IF ( B  —  X)  system. '  ’ 
Table  I  shows  the  A,,  coefficients  obtained  from  this  analysis 
and  used  in  the  spectral  fitting  procedure  to  derive  time- 
dependent  populations  from  spectra. 

For  each  vibrational  band,  the  rotational  manifold  was 
modeled  using  the  following  expression  for  the  relative  in¬ 
tensity  for  a  given  spin-rotational  transition. 


/WU(M 


e*pL~£(-v-y)/A:7"i  F 

Q,o,(T) 


( m ). 


(2) 


Terms  in  Eq.  (2)  are  defined  as  follows:  N  and  J  are  the 
upper  state  rotational  angular  momentum  and  total  angular 
momentum  excluding  nuclear  spin,  respectively;  AWand  A./ 
are  the  changes  in  these  parameters  resulting  from  emission. 
NHpf  'S"  )  closely  approximates  a  Hund’s  case  b  mole¬ 
cule,  for  which  A N  —  ±  1,  A/  =  0,  +  1  transitions  are  di¬ 
pole  allowed;  transitions  with  AAf  =  AJ  dominate,  especially 
at  larger  J.  <r  is  the  frequency  of  the  transition  in  cm  '  5  is 

the  Honl-London  or  "linestrength"  factor  appropriate  for  a 
case  b  molecule."1  The  exponential  term  is  the  Boltzmann 
factor  which  gives  the  population  in  the  upper  rovibrational 
state.  Rotational  term  values  E(NJ)  are  taken  from  Boud- 
jaadar  et  al.,'9  and  are  used  as  well  to  calculate  the  rotational 
partition  function  Qm, .  Rotational  equilibrium  is  assumed  at 
an  effective  temperature  T.  It  is  seen  in  practice  that  at  the 
total  pressures  used  in  the  LABCEDE  experiments  (i.e.,  4- 
25  Torr,  with  a  mean  collision  time  of  0. 1  fis),  rotational 
equilibrium  is  achieved  on  a  time  scale  short  compared  to  the 
time  scale  of  fluorescence  decay  ( 1-8  ms). 

,  (m)  is  the  Herman-Wallis  factor  accounting  for 
the  effect  of  centrifugal  distortion  on  the  rotational  depen¬ 
dence  of  the  line  intensity,20  where  m  =  ,V  for  R  branch  tran- 


TABLE 1.  Calculated  Einstein  coefficients  ,  for  NHl.V’I  .  i>). 


V 

1 

51.7 

2 

92.3 

3 
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sitions  and  —  (N  +  1 )  for  P  branch  transitions.  The  values 
F„„  _  ,  ( m )  were  computed  using 

=  1  +  C„„_ ,  m  +  A*-i  m2.  (3) 

in  which  the  following  values  of  the  coefficients  C„„_  ,  and 
Dn.  _  , ,  taken  from  the  semi-empirical  results  of  Chackerian 
et  al.,2n  have  been  used:  C,0  =  0.062  77,  C,,  =  0.055  11, 
C,,  =  0.048  89,  Z>10  =  0.002  251,  =  0.001  871, 

D32  =  0-  001  585.  Spectral  fits  improved  significantly  when 
the  Herman-Wallis  effect  was  thus  taken  into  account. 

In  order  to  model  the  observed  spectra,  line  intensities 
resulting  from  Eq.  (2)  were  convolved  with  the  self-apo- 
dized  instrument  lineshape  of  the  interferometer.  When  sin¬ 
gle-line  helium  lamp  radiation  is  passed  through  the  interfer¬ 
ometer.  ar.  interferogram  is  observed  which  can  be 
approximated  by  a  linear  combination  of  a  rectangularly  and 
a  triangularly  apodized  cos(.t)  function,  where  x  is  the  mir¬ 
ror  position.  The  linearity  of  the  Fourier  transformation  im¬ 
plies  that  the  instrument  lineshape  function  for  a  single-line 
source  must  be  a  linear  combination  of  a  sinc(er)  and  a 
sinc:(<7)  function.  The  extent  to  which  each  function  con¬ 
tributes  can  be  determined  by  fitting  either  the  single-line 
spectrum  or  a  molecular  emission  spectrum  such  as  NH(u), 
provided  the  lines  have  well-known  positions  and  are  ade¬ 
quately  spaced.  The  lineshape  function  depends  on  the  total 
pressure  in  the  target  chamber,  presumably  due  to  variation 
in  the  spatial  distribution  of  the  emitting  volume  with  gas 
density.  This  effect  was  taken  into  account  in  the  spectral 
fitting  that  determined  time-dependent  NH(u)  populations. 

To  determine  lineshape  parameters  for  experimental 
data,  several  NH(r)  spectra  were  first  fit  using  a  nonlinear 
least-squares  algorithm.21  In  these  fits  the  parameters 
FWHM  linewidth  A.  fraction  sine  function  a.  and  rotational 
temperature  T  were  allowed  to  vary.  Best-fit  values  of  T were 
near  room  temperature.  These  parameters  were  then  fixed  in 
computing  three  basis  sets  for  v  =  1-3,  each  basis  set  consist¬ 
ing  of  the  set  of  lines /xv  a/  convolved  with  the  instru¬ 
ment  lineshape,  then  multiplied  by  the  factor  A„.  Linear 
least-squares  fits  to  observed  spectra  using  these  basis  sets 
gives  three  multipliers,  which  are  the  relative  populations. 
Figure  2  shows  the  three  basis  functions  used  in  the  least- 
squares  fit  to  the  spectrum  in  Fig.  1.  Repeating  the  fit  at  all 
times  sampled  both  during  and  following  beam  pulse  excita¬ 
tion  of  the  gas  mixture  yields  population  decay  curves  for 
u  =  1-3,  which  can  then  serve  as  input  to  the  kinetic  analysis 
of  the  system. 

Since  both  angular  divergence  and  off-axis  entrance  of 
fluorescence  into  a  Michelson  interferometer  can  result  in 
apparent  red  shifts,  a  lineshift  parameter  <5.  proportional  to 
frequency,  was  also  introduced  into  the  nonlinear  fits.  Fitted 
values  for  5  were  <0.1%  of  the  frequency — more  typically 
0.02% — a  small  fraction  of  the  spectral  resolution,  and  too 
small  to  significantly  alfect  the  quality  of  the  fits. 

As  discussed  previously,11  spectra  not  corrected  for  the 
detection  system  response  show  equal  noise  at  all  wave-num¬ 
ber  values;  division  by  the  spectral  correction  function  re¬ 
sults  in  a  noise  level  which  increases  with  wave  number  for 
this  system.  Thus,  points  were  statistically  weighted  as  the 
inverse  of  the  spectral  correction  function  in  the  linear  fits. 
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FIG.  2.  Basis  functions  for  v  =  1-3  multiplied  by  the  population  factors 
determined  in  the  best  least-squares  fit  to  the  spectrum  in  Fig.  1 .  Parameters 
used  in  basis  functions:  lineshape  parameter  a  =  0.24,  FWHM  A  —  14.3 
cm ' rotational  temperature  T  =  300  K. 

placing  more  importance  on  data  near  the  red  end  of  the 
spectrum  where  the  detector  is  most  sensitive.  The  detection 
system  is  about  40%  Jess  sensitive  at  3400  cm  " 1  than  at  2500 
cm" '. 

The  spin-rotational  envelope  for  a  given  v  level  is  fit  well 
using  Eq.  (2),  judging  by  the  near-unity  value  of  Q  arising 
from  the  fitting  program.21  (Q  is  the  probability  that  the 
value  of  x2  resulting  from  a  spectral  fit  has  occurred  as  a 
result  of  data  uncertainty,  and  not  through  application  of  an 
incorrect  model. )  Thus  we  have  some  confidence  in  the  line- 
strength  factor/Herman-Wallis  coefficient  approach  to  the 
line  intensities.  Clearly,  the  accuracy  of  the  fitted  popula¬ 
tions  />,.  also  depends  on  the  accuracy  of  the  probabilities  A , , 
since  products  A,,  -p,  result  as  best-fit  values.  However,  the 
relaxation  rate  constants  k,,(M)  are  totally  insensitive  to  the 
absolute  values  of  the  probabilities,  and  relatively  insensitive 
to  probability  ratios.  Thus,  we  find  that  the  variation  in  rate 
parameters  calculated  when  probability  ratios  are  varied  by 
±  20%  is  well  within  the  uncertainty  arising  from  other 
sources.  The  nascent  vibrational  distribution  determined 
from  the  kinetic  analysis  is  somewhat  more  sensitive  to  the 
probability  ratios  A  r/A,f .  The  error  bars  associated  with  the 
branching  fractions  have  been  increased  to  take  this  addi¬ 
tional  uncertainty  into  account. 

IV.  KINETICS  METHODS  AND  RESULTS 

The  spectral  fits  produce  a  family  of  curves  describing 
population  evolution  of  NH(v  =1-3)  over  time.  Figure  3 
shows  these  populations  evolve  both  during  and  following 
the  beam  pulse.  NH(v)  populations  could  be  determined  at 
all  times  from  the  N2/H2  data  set.  As  mentioned  above,  in¬ 
tense  Ar*  fluorescence  prevents  observation  of  NH(v)  dur¬ 
ing  the  beam  pulse  for  the  Ar/N2/H2  data  set.  Owing  to  this 
difference,  two  different  approaches  were  taken  in  the  kinet¬ 
ic  analyses  of  the  two  data  sets;  these  approaches  are  ex¬ 
plained  in  detail  below. 

For  the  kinetic  analyses,  we  are  guided  by  the  results  of  a 
similar  analysis  for  the  OH(v)/02  system."  It  was  shown 
that  severely  correlated  and  uncertain  best-fit  parameters 


FIG.  3.  Run  consisting  of  NH(v  =  1-3)  time-dependent  populations  de¬ 
rived  from  spectral  fits  to  N./H.  data.  The  e  beam  turns  on  at  /  =  0  ms  and 
turns  off  at  /  =  0.30  ms.  While  the  beam  is  on  theNH(v  =  1-3)  populations 
are  seen  to  increase:  decays  following  beam  termination  are  caused  mainly 
by  colltsional  relaxation.  Also  shown  is  a  kinetic  fit  (solid  curves)  to  numer¬ 
ous  runs  of  this  type,  which  determines  k,a  (NO  and  relative  creation  rates 
for  v  =  1-3.  Pressures  used:  15.2  Torr  N\.  0.20  Torr  H,. 

result  when  more  than  one  kinetic  parameter  is  introduced 
for  each  v  level  in  the  fit.  Thus,  in  practice  it  is  not  feasible  to 
accurately  determine  kinetic  parameters  describing  all  possi¬ 
ble  decay  channels,  including  multi-quantum  decay  and  re¬ 
active  loss.  Studies  of  other  systems  such  as  OH(v  =  1-3)," 
NO(v=  l-7),J2andCO(v  =  1-1 6)"' have  shown  that  a  sin¬ 
gle-quantum  mechanism  can  be  used  to  successfully  model 
vibrational  relaxation  in  diatomic  molecules;  it  is  the  one  we 
use  here. 

No  exothermic  pathways  exist  for  the  reaction  of 
NH{v<;3)  with  the  bath  gases  Ar  and  N;.  This  observation 
is  consistent  with  the  measured  upper  limit  of  3(  —  19) 
cm ’/s  for  the  rate  constant  describing  reaction  between  NH 
and  N2.2'  In  the  case  ofH,,  however,  the  reaction 

NH  +  H--NH.  +  H. 

AH"  =  0.61  ±0.06  eV  (Refs.  24.25)  (4) 

becomes  exothermic  for  NH(v>2).  Thus,  while  reaction 
(4)  does  not  occur  for  vibrational  ground  state  NH 
[k<H  —  17)  cmVs  (Ref.  26)  J,  it  could  play  a  significani 
role  for  higher  v  levels. 

The  electron  beam  bombardment  of  N:  and  H;  also 
gives  rise  to  species  which  are  not  directly  monitored  in  our 
experiment,  notably  NH(v  =  0)  and  N  and  H  atoms.  All  of 
these  species  react  rapidly  with  NH: 

NH  +  NH  — products,  8.5(  —  1 1  )cm3/s(Ref.  27). 

(5) 

NH +  N-R.  +  H.  1.8(  —  1 1  )cmVs(Ref.  3). 

5.0(  —  ll)cmVs(Ref.  28).  (6) 

NH  +  H-Hj  +  N,  3.0(  —  12)cmVs(Ref.3), 

1.66(  —  10)cmVs(Ref.  28),  (7) 

NH(r)  +  M-NH(n-l)  +  M,  M  =  H.N.  (8) 

Spencer  and  Glass2'1  measured  the  rate  constant  for  deacti¬ 
vation  of  OH(u  =  1.2)  by  H  and  found  it  to  be  gas  kinetic; 
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the  relaxation  of  NH(y)  by  H  atoms  may  also  be  fast.  Pre¬ 
vious  work  on  the  NH  system5  has  shown  that  ionic  path¬ 
ways  lead  predominantly  to  the  production  of  H  atoms. 
Also,  as  demonstrated  by  our  modeling  studies  (see  below), 
the  populations  of  the  beam-created  species  are  complicated 
functions  of  pulse  width,  beam  current,  duty  cycle,  and  gas 
pressure.  In  this  work  we  have  used  conditions  which  should 
minimize  build-up  of  these  quenchers  in  the  target  chamber. 

Ideally,  the  data  allow  relaxation  rates  to  be  determined 
independently  for  all  three  observed  u  levels.  Population  in 
NH (o  =3),  however,  was  present  in  relatively  small  quan¬ 
tity  compared  with  NH(d  =  1,2),  with  poorer  S/N  (Fig. 
3 ) .  Because  of  the  small  signal  for  v  =  3.  relaxation  rates  for 
v  =  3  could  not  be  determined  from  the  data  with  acceptable 
precision.  On  the  other  hand,  the  effect  of  feed  from  u  =  3 
into  lower  v  levels  is  significant  and  must  be  taken  into  ac¬ 
count  in  the  empirical  model.  In  order  to  fit  the  data  with 
fewer  parameters,  a  power  law  dependence 
=  u"k,.  „  |  (M)  was  assumed  for  each  relaxer  gas  M.  Such  a 
dependence  has  been  found  to  hold  for  numerous  systems.'0 
and  represents  a  minimal  constraint  on  the  empirical  model. 
The  adjustable  parameters  k„  „  ,  and  a  are  largely  deter¬ 
mined  by  the  u  =1.2  levels,  which  have  greater  S/N.  Rate 
constants  fort,1  =  1.2  are  thus  determined  with  greater  preci¬ 
sion  than  those  for  v  =  3. 

Data  were  fit  to  both  NH  growth  and  decay  parameters 


using  the  formulae: 

^y111  =  T-.ftO  -*-.,[NH(l)l 
a( 

+  *h.jINH(2)], 

(9) 

d  [  NH(  2)  ]  =x^:p{l)  _fir_,[NH(2)] 
at 

+  j  [NH(3)  ], 

(10) 

d  [ NH ( 3 ) ]  =  _R"_  JNH(3)]. 

at 

(ID 

The  first  term  in  Eqs.  ( 9 )— ( II)  describes  the  effects  of 
chemical  feed  from  the  N(:Z>)  +  H-  reaction;  =  i-i  are 
coefficients  giving  the  relative  feed  rate  into  NH(v  =  1-3), 
and  Pit)  is  a  function  describing  the  time  dependence  of  the 
precursor  concentration;  it  is  discussed  in  more  detail  below. 
Because  the  function  Pit)  includes  a  multiplicative  factor 
(see  below),  the  four  parameters  %vm , _3  and  Pit )  include  a 
redundancy;  thus,  in  the  fits  ,  is  set  to  unity  and  „  2  J 
are  solved  as  fractions  of  the  u  =  1  feed  rate.  Relative  cre¬ 
ation  rates  for  NO(t>  =  1-7)  have  been  analyzed  in  a  similar 
fashion." 

The  second  and  third  terms  represent  decay  from  and 
feed  into,  respectively,  a  given  level  v.  Single-quantum  decay 
rate  parameters  R„  are  defined  by 

R„  =  ,  (M)  [M]}  -1-  R,,„,  (12) 

where  the  sum  is  over  all  of  the  bath  gases  (Ar,  N:,  H, ) 
present  in  the  reaction  chamber.  Here,  Rr0  is  residual 
NH(v  =  1-3)  loss,  due  largely  to  quenching  by  beam-creat¬ 
ed  species,  which  we  assume  to  be  a  constant  over  the  gas 
pressures  used  in  these  experiments.  This  assumption,  while 
only  an  approximation,  is  supported  by  the  results  of  both 


the  empirical  fits  and  the  kinetic  modeling  (see  below).  Sin¬ 
gle-quantum  radiative  rates  Av_„_  ,  contribute  only  to  the 
quantities  R„n,  and  double-quantum  radiative  rates  , 
are  negligible  compared  to  other  terms. 

The  time  dependence  of  the  N  ( 2D)  concentration  is  giv¬ 
en  by  Pit),  mentioned  above.  Assuming  that  the  rate  of  for¬ 
mation  of  N(_D)  is  constant  during  the  e  beam  and  zero 
after  beam  turn-off,  Pit)  is  given  by 

Pit)  =  C[1  —  exp( -R,t)]  (13) 

while  the  e  beam  is  on.  and 

Pit)  =  C  exp[  —/?,(/  —  **)]  (14) 

while  the  e  beam  is  off,  where  the  beam  turns  on  at  t  =  0  and 
turns  off  at  t  =  th,  and  C  is  a  constant.  The  precursor  decay 
rate  Rn  is  given  by 

/?r  =  {S^(M)[Ml}  +  /?,.0,  (15) 

where  the  sum  is  over  the  bath  gases  Ar,  N;,  and  H:,  and  Rn0 
describes  residual  N (:Z>)  loss  which  we  assume  to  be  a  con¬ 
stant  with  respect  to  gas  pressure.  The  following  values  for 
kn(M)  were  used  (units  of  cmVs):  ^(H-)  =  2.3(  —  12);7 
kp <N:)  =  1.5(  -  14>;,J jfc„(Ar)  =  1.0(  -  16)." 

The  rates  7?..  and  the  precursor  population  Pit)  are 
functions  of  gas  partial  pressures,  and  thus  depend  on  the 
particular  data  run  being  examined.  However,  the  rate  con¬ 
stants.  k,.  _  ,  ( M)  and  a(M)  (M  =  Ar.N:,H;).  the  creation 
rate  multipliers  =  ,_3 ,  and  the  residual  quenching  rates 
R,.{)  and  Rp0  are  independent  of  individual  run  conditions. 
Therefore,  a  “global”  fitting  technique  was  used  in  which  the 
time-dependent  populations  were  fit  simultaneously  for 
these  parameters.  In  addition,  an  adjustable  gain  parameter 
C  is  required  for  each  run.  Such  a  fit  was  necessitated  by  our 
use  of  k„  „  ,  and  a  as  decay  rate  parameters,  since  a  fit  to  an 
individual  run  cannot  determine  these  two  parameters.  For  a 
given  run  the  coupled  linear  differential  equations  were 
solved  numerically  by  estimating  initial  values  for  the  rate 
parameters,  then  integrating  the  equations  from  the  time  of 
e-beam  onset  until  the  time  of  the  last  data  point  to  produce 
model  curves  for  v  =  1-3.  A  nonlinear  least-squares  fitting 
routine  then  iteratively  adjusted  the  parameters  to  produce 
the  best  least-squares  fit  to  the  experimental  data.  The  pres¬ 
ent  approach  takes  the  place  of  the  two-step  analysis  used  in 
our  previous  study  of  OH(u)  relaxation  by  0:.H  in  which 
data  were  first  fit  on  a  run-by-run  basis,  and  then  rate  con¬ 
stants  determined  as  the  slope  of  the  resultant  rates  plotted 
against  O-  partial  pressure. 

Because  the  N-/H,  and  Ar/N-/H-  data  sets  were  taken 
under  different  conditions,  the  two  data  sets  were  fit  individ¬ 
ually  to  appropriate  subsets  of  the  global  parameters  de¬ 
scribed  above.  Excellent  fits  resulted  when  this  procedure 
was  followed.  The  two  fits  are  described  below. 

The  N,/H;  data  set  was  first  fit  using  ten  global  param¬ 
eters  ,(N,),  a(N:),  ,  (HO,  a(H-). 

IV-  i.i  <  I-,.,,,  and  Because  the  resulting  values  of 

the  H-  kinetic  parameters  k„  were  not  different  from  zero 
within  systematic  error,  the  fc„(H:)  were  set  equal  to  zero 
and  the  fit  repeated.  A  total  of  16680  data  points  ( 139  runs, 
each  consisting  of  40  population  values  for  each  of  three 
vibrational  levels  at  25  fi s  intervals)  were  fit  using  the  rc- 
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maining  eight  global  parameters.  A  good  fit  resulted  from 
this  model,  as  evidenced  in  the  4.3%  rms  deviation  of  the  fit 
relative  to  the  largest  data  point  in  a  given  run,  which  is 
within  the  estimated  5%  uncertainty  in  the  average  data 
point.  For  N2.  kv_  ,  =  1.5(  -  14)  cmVs  and  a  =  1.4  were 
determined  from  the  fit  to  the  N2/H2  data  set.  The  normal¬ 
ized  relative  formation  rates  of  NH(e  =  1-3 )  were  found  to 
be  0.36,  0.35,  and  0.29,  respectively.  Best-fit  relative  forma¬ 
tion  rates  can  be  equated,  in  turn,  to  the  nascent  vibrational 
distribution  in  NH(n=  1-3)  resulting  from  N(:Z>)  +  H:. 
These  measured  branching  fractions  are  given  in  the  first 
column  of  Table  II.  Figure  3  shows  the  results  of  this  fit 
applied  to  a  particular  run. 

In  the  second  data  set,  where  Ar/N,/H:  mixtures  were 
used,  data  were  taken  after  termination  of  the  e  beam,  and 
therefore  the  NH  ( v )  relative  creation  rates  were  not  used  as 
adjustable  parameters.  In  the  absence  of  literature  values  for 
the  creation  rates,  the  above  rates  were  used  as  fixed  param¬ 
eters  in  the  fit.  Best-fit  values  for  A,a  ( Ar.N;)  were  found  to 
vary  by  less  than  5%  upon  50%  variation  of  the  relative 
creation  rates.  The  N(:D)  decay  parameter  Rp,t  could  not 
be  determined  and  was  set  equal  to  5x  10'  s'  ',  the  value 
obtained  from  the  fit  to  the  N2/H2  data.  The  fit  to  the 
Ar/N2/H2  data  also  gave  poorly  determined  H;  kinetic  pa¬ 
rameters;  thus,  A,.  ( H2 )  were  set  equal  to  zero  and  the  fit 
repeated.  A  total  of  73308  data  points  (222  runs,  consisting 
of  an  average  of  1 10  populations  for  each  of  three  vibrational 
levels  at  10 /xs  intervals)  were  fit  using  seven  global  param¬ 
eters  [ k,. ,  ,  (Ar).  a(Ar).  k, ,  (N;).  o(N:),  R, . 

The  fit  resulted  in  an  rms  deviation  of  3.9%,  similar  to  the 
rms  deviation  resulting  from  the  analysis  of  the  N:/H:  data. 
Values  of  k ,  (Ar)  =  2(  —  15)  enr/s.  a(Ar)  =  1.2.  and 
A,._  ,  (N2)  =  1.0(  —  14)  em'/s.  cr(N;)  =  l.°  were  deter¬ 
mined  from  the  fit  to  the  Ar/N-/H-  data  set.  Figure  4  shows 
the  fit  superimposed  on  data  points  for  a  particular  run. 

Final  values  for  A,.„  (N;)  were  chosen  to  be  the  aver¬ 

age  of  the  results  from  the  N,/H-  and  Ar/N2/H;  studies, 
since  the  two  data  sets  were  composed  of  a  similar  number  of 
runs,  and  because  the  fits  were  of  comparable  quality.  While 
the  values  of  the  A,  (N\)  resulting  from  the  two  fits  differ  by 
more  than  their  combined  statistical  uncertainties,  the  dif¬ 
ferences  arc  within  the  estimated  systematic  error  in  the  ex¬ 
periment.  Neither  data  set  produced  values  for  H;  kinetic 
parameters.  Instead,  upper  limits  for  A,  (H.)  were  set  at  the 
maximum  values  which  produced  acceptable  fits.  Using  this 
method.  A, (H2)<0.5,  1.0,  and  1 .5 (  —  12)  em'/s,  re¬ 
spectively,  where  a  has  been  assumed  to  be  1  0.  All  of  the 
relaxation  kinetics  results  are  collected  in  Table  III. 


Values  of  R,  0  =  500-5000  s"  1  were  obtained  in  the  fits 
to  the  N2/H2  and  Ar/N,/H:  data  sets.  These  values  are  sig¬ 
nificantly  larger  than  can  be  accounted  for  by  either  radia¬ 
tive  loss  or  diffusion  out  of  the  field-of-view  (10  s"  1  at  the 
lowest  pressures  used).  Instead,  these  values  probably  re¬ 
flect  quenching  by  beam-created  species.  NH(u)  is 
quenched  to  some  degree  by  itself  and  by  H  and  N  [Eqs. 
( 5 )— ( 8)  ] .  The  N  (2D)  residual  decay  parameter 
Rr0  =  5000  s  ~ 1  determined  in  the  fit  to  the  N2/H2  data  set 
is  also  consistent  with  efficient  quenching  of  N(:Z3)  by 
beam-created  species,  which  could  include  electrons,  ions, 
and  NH. 

To  investigate  the  role  of  beam-created  quenchers,  a 
number  of  runs  were  taken  under  various  conditions  of  e- 
bcam  pulse  width  (250-2750 g/s)  and  duty  cycle  period  (4- 
30  ms),  with  constant  N-  and  H;  partial  pressures.  The 
NH(r>0)  relaxation  rate  increases  linearly  with  increased 
pulse  width,  apart  from  an  offset  due  to  N-  and  H-  quench¬ 
ing.  This  result  suggests  a  proportionality  relationship 
between  pulse  width  and  the  steady-state  concentration  of 
beam-created  qucnuiers.  The  slope  implies  beam-induced 
quenching  qualitatively  consistent  with  the  residual  NH  ( r ) 
loss  rates  determined  in  the  empirical  fits.  Assuming  a  gas- 
kinetic  quenching  rate  constant  of  10'  em'/s,  a  residual 
loss  rate  of  1000  s ' 1  implies  a  steady-state  quencher  concen¬ 
tration  of  101'  cm  ‘  '.  in  qualitative  accord  with  expectation. 
Larger,  but  still  reasonable,  quencher  concentrations  are  cal¬ 
culated  if  one  assumes  smaller  quenching  rate  constants. 
The  above  study  determined  that  the  conditions  which  mini¬ 
mize  beam-created  quencher  activity  are  pulse  width  CbOO 
fj s  and  duty  cycle  period >7.5  ms.  Within  this  regime,  the 
empirical  model  provides  a  satisfactory  description  of 
NH(c)  decay  due  to  residual  quenching.  All  data  used  to 
determine  rate  constants  were  collected  under  toese  con¬ 
straints. 

V.  MODELING  OF  REACTION  CHEMISTRY 

In  order  to  test  hypotheses  concerning  formation  mech¬ 
anisms  and  reactivity  of  beam-created  species,  we  have  mod¬ 
eled  the  chemistry  of  irradiated  N;/H;  mixtures  using  a  de¬ 
tailed  chemical  model.  The  calculations  use  subroutines 
from  the  CHEMKIN  chemical  kinetics  software  package.'"1 
together  with  estimates  of  the  initial  yields  of  the  different 
reactive  species,  both  neutral  and  ionic,  created  by  the  inter¬ 
action  of  the  electron  beam  with  the  gas  mixture.  The  initial 
yields  were  estimated  using  the  Electron  Energy  Loss 
Code'*  to  determine  the  effects  of  the  deposited  energy.  The 


TABLE  II  Naccent  product  distribution  for  N(:£>)  —  H-  —  NH(A-'L  .i"'3)  -  H 


r 

Measured 

Branching  fraction 
Prior' 

Measured  -1-  Surpnsaf' 

0 

0.53 

0  15 

t 

0.36  *  0.10 

0.30 

0.3)  -  0.09 

2 

0.35  i  0.10 

O.H 

0.30  -  0.09 

3 

0.29  ±0.10 

0.03 

0.24  -  ()  0« 

•Equation  (18) 

•Branching  fraction  for  v  =  0 determined  from  Fig  6(b). 
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FIG.  4.  Run  consisting  of  NH(v  =  1-3)  time-dependent  populations  ob¬ 
tained  from  Ar/N./H,  data.  The  e-beam  turns  on  at  t  =  0  ms  and  turns  ofT 
at  i  =  0.60  ms.  Also  shown  is  a  kinetic  fit  (solid  curves)  to  numerous  runs  of 
this  type,  which  determines  k.a  ( Ar.N.).  Pressures  used:  14.3  Torr  Ar.  2.') 
Torr  N.,  0.08  Torr  H,. 


interaction  volume  was  estimated  from  the  relation  for  the 
cone  angle1'’  in  pure  N:,  while  providing  for  the  small  effect 
of  added  H:  using  scaling  based  on  the  ionization  cross  sec¬ 
tion.  Model  calculations  typically  included  100  reactions. 17 

The  kinetic  model  allows  species  to  evolve  in  time  over 
multiple  duty  cycles  assuming  specified  parameters,  e  g.,  gas 
partial  pressures,  duty  cycle  period,  and  e-beam  pulse  width 
and  current.  The  most  important  beam-created  species  are  H 
atoms  and  N  atoms,  both  of  which  reach  steady-state  con¬ 
centrations  over  several  e-beam  duty  cycles.  The  maximum 
concentrations  of  the  H  atoms  are  2-10  times  the  maximum 
N  atom  concentrations.  For  example,  in  a  mixture  of  7.5 
Torr  N.  and  0.05  Torr  H-,  with  a  300  /xs  pulse  width,  one 
determines  [H]  =4(13)  cm-3,  [N]  =2(13)  cm'3,  and 
[NH]  =3(10)  cm'3.  These  concentrations,  when  multi¬ 
plied  by  the  rate  constants  for  reactions  ( 5 )— ( 7 ) ,  give  480- 
7640  s'  1  for  the  contribution  of  H,  N,  and  NH  to  NH(u) 
quenching,  depending  on  the  choice  of  the  rate  constants. 
NH(u)  quenching  by  NH  is  negligible  compared  to  quench¬ 
ing  by  H  and  N  atoms.  The  480  s  "  1  value  is  a  lower  limit, 
since  we  have  neglected  any  contribution  from  reaction  (8). 
Moreover,  the  lower  value  for  the  rate  constant  for  reaction 
(7)  includes  the  effects  of  an  activation  energy  of  1000  K,  a 
barrier  that  may  decrease  upon  vibrational  excitation  of  the 
NH  reactant.  Thus,  estimated  rates  for  the  quenching  of 
NH(u)  by  H  and  N  are  consistent  with  fitted  values  Ro0. 

The  kinetic  model  generates  both  relaxation  rate  con- 


TABLE  III.  Rate  constants  for  NH(.T'Z',  i>)  M  — NH(u  —  I )  4-  M. 


u 

Rate  constant  ( 10  14  cm'  s 
N,  Ar 

') 

H,4 

1 

1.2  ±0.3 

0.2  ±0.1 

<  50. 

2 

3.3  ±  1.5 

0.5  ±0.2 

<100. 

3 

7.5 ->-2.5 

0.8  ±0.3 

<150. 

'  May  include  reactive  quenching  via  reaction  (4)  for  NH(r  =  2.3). 


stants  and  ratios  of  NH(u)  creation  rates.  Figure  5  shows 
values  for  NH(u  =1-3)  time-dependent  populations  calcu¬ 
lated  by  the  kinetic  model,  compared  with  experiment.  [Ab¬ 
solute  NH(u)  concentrations,  determined  from  the  model, 
are  consistent  with  expectation.  ]  The  model  predicts  that 
among  the  potential  sources  for  NH(u),  N(2Z>)  plays  a 
much  more  significant  role  than  other  potential  sources, 
such  as  ionic  species  and  metastable  N2(/1 32ut ).  With  rea¬ 
sonable  choices  for  the  kinetic  and  initial  yield  parameters, 
the  kinetic  model  predicts  N,  and  H,  relaxation  rate  con¬ 
stants  and  a  nascent  distribution  of  NH(u)  that  are  within 
the  limits  quoted  in  Tables  II  and  III.  These  observations 
provide  further  evidence  that  the  observed  NH(u)  result 
from  N(:Z))  4-  H2. 

It  is  important  to  draw  a  distinction  between  the  meth¬ 
odology  of  the  kinetic  model  and  that  of  the  empirical  fitting 
process.  The  empirical  fit  does  not  include  a  detailed  set  of 
reactions,  and  thus  does  not  address  the  possible  role  of 
beam-created  quenchers,  other  than  yielding  effective  values 
of  R„0  and  R n0.  The  kinetic  model,  on  the  other  hand,  in¬ 
cludes  numerous  reaction  pathways.  While  a  number  of 
these  reactions  have  poorly  determined  rate  constants,  vari¬ 
ation  of  these  rate  constants  within  reasonable  limits  does 
not  strongly  affect  predictions  made  for  rates  of  NH(e)  for¬ 
mation  and  decay.  As  a  result,  the  model  remains  in  qualita¬ 
tive  agreement  with  the  data  for  a  wide  range  of  chosen  input 
parameters.  Model  predictions  did  depend,  however,  on  the 
inclusion  of  pathways  describing  quenching  of  NH(y)  by 
beam-created  species,  especially  H  atoms.  Thus,  the  model 
gives  insight  into  the  role  of  the  beam-created  N  and  H 
atoms  on  the  population  of  NH(u>0).  It  also  provides  a 


FIG.  5.  Semiloe  plot  showing  kinetic  model  for  NH(v  =  1-3)  time  evolu¬ 
tion  superimposed  on  (normalized)  data.  The  e  beam  turns  on  at  l  =  0  ms 
and  turns  ofTat  r  =■  0.30  ms.  Absolute  concentrations  for  NH(u)  were  de¬ 
termined  from  :h^  model  and  not  from  experiment,  which  provided  only 
relative  populations. 
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valuable  check  on  possible  problems  in  the  empirical  fitting 
process,  which,  if  present,  would  give  rise  to  poorer  agree¬ 
ment. 


VI.  STATISTICAL  ANALYSIS  OF  PRODUCT 
VIBRATIONAL  DISTRIBUTION 

Because  we  detect  product  NH(u)  in  emission,  we  can¬ 
not  observe  formation  of  NH(u  =  0).  However,  the  nascent 
yield  for  NH(u  =  0)  can  be  inferred  by  comparing  our  re¬ 
sults  for  NH(i/=  1-3)  with  a  distribution  obtained  by  as¬ 
suming  a  statistical  partitioning  of  the  total  energy  available 
to  the  NH  ( v)  product.  The  statistical  model  assumes  that  all 
product  quantum  states  are  equally  probable,  given  the  sole 
constraint  of  conservation  of  total  energy.  This  surprisal 
analysis  method  has  been  discussed  in  numerous  reports.”1 

To  make  a  prediction  for  the  relative  amount  of 
NH(i>  =  0),  the  surprisal  S(/. )  is  plotted  as  a  function  of  the 
reduced  energy  /  =  £„/<£r),  where  £,  is  the  energy  of 
vibrational  state  NH(t>)  relative  to  NH ( r  =  0),and  < ET )  is 
the  average  total  energy  available  to  the  NH(t')  product. 
Thus,/  is  the  fraction  of  the  available  product  energy  pres¬ 
ent  in  vibration.  £,  was  calculated  for  v  =  0-3  using  the 
band  origins  determined  in  the  spectroscopic  study  by  Boud- 
jaadar  et  a/.'1’  ( ET )  is  estimated  as 

(Er)  =AE"  +  Ea  +lkT+kT,  (16) 

where  A£"  and  £„  are  the  reaction  energy  and  activation 
energy,  respectively.  The  sum  £„  4-  (3/2) A:  7"  is  the  average 
reactive  collision  energy,  and  the  last  term  is  the  average 
rotational  energy.’1’  (£r>  is  uncertain  to  ±  kT. 

To  determine  A£",  we  use  the  energy  cycle 

N(:Z»  +  H:  — N(4S)  4-  H--N(4S)  4-  2H 

-NHl’S'l-f-H.  (17) 

A £°(eV)  =  -  2.3846(  1 )  4-  4.4781  (2)  -  3.37(3) 

=  -  1.28(3)  eV, 

with  standard  values  for  DntH-)1'  and  £(N(JS) 
—  N  i'D)). 4,1  and  the  recently  calculated  value  for 
Z>n(NH(  ’Z)).4'  The  activation  energy  £„  in  Eq.  (16)  is  tak¬ 
en  to  be  1  kcal/mole,  on  the  basis  of  measured  activation 
energies  for  hydrogen  atom  transfer  reactions  with  similar 
AH".4:  Thus  (£r>  =1.39  ±  0.04  eV  is  the  average  total 
energy  available  to  the  nascent  NH(i>)  product  from  N(:Z>) 
4-  H;. 

To  perform  a  surprisal  analysis,  a  “prior”  function 
/’"(/. )  must  be  defined  which  describes  the  prior  expectation 
for  the  NH(  r)  product  distribution  assuming  statistical  dis¬ 
tribution  of  the  available  energy.  It  can  be  shown”  that  for  a 
three-body  reaction  such  as  N  (:D)  4-  H-  the  prior  function 
can  be  approximated  by 

£"(/ ;.)  =  (1  -/.)'/:/£  (1  -/.)’",  (W 

«■ 

where  the  sum  in  the  denominator  is  over  all  populated  prod¬ 
uct  t'  states,  i.e..  v  =  0-3.  Similar  functional  dependencies 
result  even  when  constraints  arising  from  the  conservation 
of  angular  momentum  are  taken  into  account.”1  The  prior 
values  calculated  using  Eq.  (18)  are  plotted  along  with  the 


measured  values  for  v  =  1-3  in  Fig.  6(a). 

Figure  6(b)  shows  a  plot  of  the  surprisal  function 

£</.)=  -\n[P{J„)/P°  {/„)],  ‘  (19) 

vs./,,  where  P(f, .)  is  the  observed  vibrational  distribution. 
The  relationship  is  seen  to  be  linear,  consistent  with  a  single 
reaction  mechanism  giving  rise  to  observed  NH  ( v) .  The  sur¬ 
prisal  plot  yields  NH(0):NH(1 )  =  0.47,  suggesting  that 
NH(t)>0)  levels  are  populated  more  quickly  than  the 
ground  vibrational  state.  The  middle  column  of  Table  II 
shows  the  prior  distribution  of  nascent  NH(u),  while  the 
right-hand  column  gives  the  distribution  inferred  from  mea¬ 
surement  plus  the  surprisal  analysis. 

The  negative  slope  (  —  3.9)  in  Fig.  6(b)  is  indicative  of 
a  product  vibrational  distribution  which  is  more  excited 
than  that  predicted  by  the  statistical  model,  i.e.,  the  mea¬ 
sured  average  product  level  (v)mr„  is  larger  than  the  prior 
prediction  (r)pr,or.  Here,  (v)  is  determined  from 

<„>  =  y  iTi„  (20) 

i- =  O 

where  \'r  is  the  normalized  branching  fraction.  With  the 
data  in  Table  II,  Eq.  (20)  can  be  used  to  calculate 
=  1.6  and  (t*)r„„r  =  0.7.  In  addition,  the  fraction  of 
the  average  total  product  energy  in  the  NH(i>)  vibration  is 
given  by 


(/,)  = 


<£,.) 

<£r> 


1 

( Et ) 


I  ErXr. 


(21) 


where  (£,.)  is  the  average  vibrational  energy.  With  Eq.  (21) 
and  Table  II.  one  determines  (/. )  =  0.44,  substantially  larg¬ 
er  than  the  prior  value  <£,.  )pnnr  =0.18. 


0.00  0.25  0.50  0.75  1.00 


U 

FIG.  fi.  ( a  1  Calculated  and  observed  nascent  vibrational  distribution  as  a 
function  of  the  reduced  energy  parameter/,  the  fraction  of  the  available 
produci  energy  in  the  NH  vibration  Circles  are  data  for  \  =  l-.t  and  trian¬ 
gles  arc  calculated  prior  values  for  ra  0  —  3  from  Eq  (18).  The  circle  ai 
f,  =  0  is  the  prediction  for  s  =  0  made  by  the  surprisal  analysis,  (b)  Sur¬ 
prisal  function  —  In  (P/P")  plotted  against/  The  linear  relationship  is 
consistent  with  a  single  operative  mechanism  for  NHli1)  production;  the 
negative  slope  reflects  a  product  distribution  which  is  more  highly  excited 
than  the  pnot  expectation.  The  v  intercept  gives  rise  to  the  prediction  for 
r  =  0  shown  in  cig.  6(  a ) . 
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II.  DISCUSSION 

Much  attention  in  the  field  of  vibrational  energy  trans- 
:r  has  focused  on  the  relationship  between  the  single-quan- 
im  relaxation  rate  constant  ,  and  v.  Qualitatively  dif- 
:rent  behavior  has  been  observed  for  the  vibrational 
‘taxation  of  diatomic  and  polyatomic  molecules.  In  gen- 
ral,  diatomics  relax  according  to  kv  =  uakl>  _ , ,  where  a  >  1, 
rhile  polyatomics  relax  according  to  k„  =  v  kum  ,  ,30  This 
ehavior  has  been  seen  for  relaxation  of  the  open  shell  spe- 
ies  NO(u  =  1-7),22  where  a  i/1 5  dependence  was  observed, 
s  well  as  for  relaxation  of  such  closed-shell  species  as  HF 
nd  CO.30  For  relaxation  of  NH(u  =  1-3)  by  N2  and  Ar  we 
etermine  a  values  of  1.7  ±  0.3  and  1.2  +  0.3,  respectively, 
idicating  behavior  intermediate  between  linear  and  qua- 
ratic  dependence  on  v.  Measurements  o(ku>}  will  require  a 
afferent  source  for  NH(u). 

Relaxation  of  NH(u  =  1-3)  by  both  N:  and  Ar  is  ex- 
remely  inefficient;  on  average,  about  104  collisions  are  re¬ 
tired  for  relaxation  by  N2,  and  5  X  104  collisions  for  relaxa- 
ion  by  Ar.  An  efficient  V-V  transfer  mechanism  is  not 
vailable  for  relaxing  NH(u  =  1-3)  by  collision  with  either 
4,  or  Ar  [  v(N,)  =  2330  cm  '  ',  v(NH(w  =  1-3)) 
=  3570-3230  cm'1].  Thus,  a  V-R  and/or  V-T  mechanism 
s  reflected  in  the  very  low  observed  efficiencies,  especially  in 
he  case  of  Ar. 

The  nascent  distribution  parameter  (/].)  can  be  com¬ 
pared  with  those  arising  from  the  reactions  of  F  atoms  with 
ivdrides  RH.  for  which  there  are  abundant  data.  In  general, 
eactions  F  4-  HR  — FH  4-  R  channel  a  percentage  of  the 
ivailable  energy  into  the  product  (HF)  vibration  which  is 
:omparable  to  that  of  N(:Z>)  4-  H2.  For  instance,  the  reac- 
ions  F  +  HR  with  group  VI  hydrides  HR  =  H.O,  CH.OH. 
ind  H:S  yield  values  (/„)  of  0.41, 0.44,  and  0.45,  respective- 
y.42  In  addition,  values  (/„ )  do  not  seem  to  depend  strongly 
)n  the  nature  of  R;  it  has  been  shown  that  for  reactions  of  this 
:ype  (A  4-  HB  — AH  4-  B.  with  A  and  B  both  much  heavier 
han  H),  a  relatively  small  percentage  of  the  available  pred¬ 
ict  energy  is  channeled  into  the  vibrational  energy  of  the  R 
fragment. 

The  reaction  N(2Z))  4-  H-  can  also  be  compared  with 
reaction  (22),  which  involves  an  analogous  reaction  of 
N(2Z))  with  a  homonuclear  diatomic  molecule: 

N(2Z?)  4-  02  — NO(i/<  18)  4-  0(}P),  AEn  =  -3.76eV, 

(22a) 

-NO(u<8)  4-OC0),  AE°  -  —  1.79  eV, 

(22b) 

Reaction  (22)  is  complicated  by  competing  product  chan¬ 
nels  giving  the  ground  and  first  excited  states  of  O.  respec¬ 
tively.  A  recent  analysis,43  however,  suggests  that  reaction 
(22b)  is  the  more  important  channel  for  NO(u)  formation. 
Observed  peaking  toward  high  u  levels  in  the  product  implies 
a  direct  abstraction  mechanism  with  a  repulsive  approach 
configuration  and  a  short-lived  collision  complex.  The  reac¬ 
tion  of  higher  energy  N(  2P)  with  02,  on  the  other  hand,  may 


occur  via  a  long-lived  complex  with  complete  randomization 
of  the  reaction  energy.43  In  this  case  large  values  of  J  in 
initially  formed  product  states  NO(uJ)  are  consistent. with 
an  attractive  approach  configuration. 

Further  comparisons  with  N  ( 2D)  4-  H2  can  be  made  by 
considering  reactions  (23 )— (25 ),  all  of  which  nominally  in¬ 
volve  hydrogen  atom  transfer  between  a  second-row  atom 
and  H2: 

N(2Z?)  4-  H;-NH(k3)  4-  H.  AE°  =  —  1.28  eV, 

(23) 

O CD)  +  H-  — OH(u<4)  4-  H.  AE°  =  -  1.88  eV, 

(24) 

FC-P)  4-H2~FH(o<3)  4-H,  AE°  =  -1.38  eV. 

(25) 

F  4-  H:  —  FH  4-  H  is  one  of  the  most  closely  studied  atom 
transfer  processes.  The  reaction  produces  a  vibrational  dis¬ 
tribution  which  is  strongly  peaked  at  u  =  2.44  Most  of  the 
available  product  energy  is  channeled  into  FH  vibration: 
(/],)  =  0.67. 44  These  observations,  along  with  trajectory 
studies,  suggest  that  F  4-  H2  proceeds  as  a  direct  atom  ab¬ 
straction.43  Reactions  F  4-  HR  — FH  4-  R  mentioned  above 
probably  involve  the  same  mechanism. 

The  product  level  distribution  for  O('D) 
4-  H:  —  OH(d)  4-  H  has  also  been  discussed  in  many  re¬ 
ports46  with  some  remaining  disagreement.47  Studies  indi¬ 
cate  a  broad,  fairly  constant  distribution,  with  measured  val¬ 
ues  of  (/„>  =  0.334S  and  (/„>  =0.3949  Analogous  to  the 
reaction  N(:/>)  4-  02  —  NO(u,J)  4-  O,  mentioned  above, 
0('£>)  4-  H2  generates  product  states  OH(y,/)  with  large 
values  of  J.  Molecular  symmetry  considerations  show  that 
theO('Z))  4-  H:  reaction  can  occur  on  two  different  surfaces 
to  give  product  species,  one  surface  involving  direct  hydro¬ 
gen  atom  transfer,  and  one  surface  a  metastable  intermedi¬ 
ate.6  The  second  mechanism,  using  0(  '£>)  insertion  to  form 
a  short-lived  H-O-H  intermediate  which  then  dissociates,  is 
thought  to  be  the  predominant  channel.47 

For  N(:Z?)  4-  H2  —  NH(u)  4-  H,  the  observed  linear 
surprisal  plot  is  consistent  with  a  single  formation  mecha¬ 
nism.  Unlike  reported  product  i>  distributions  of  the 
O('D)  4-  H,  reaction,  the  NH(u)  product  distribution  (Ta¬ 
ble  II)  disfavors  u  =  0  product.  Also,  the  average  NH(t>) 
product  energy  (/„)  =  0.44  is  similar  to  that  of  the  HF(u) 
product  of  F  4-  HR  reactions,  discussed  above.  These  obser¬ 
vations  suggest  a  direct  abstraction  mechanism,  as  opposed 
to  a  mechanism  that  involves  an  H-N-H  intermediate. 
Clearly,  however,  further  study  is  needed  in  order  to  eluci¬ 
date  the  mechanism  for  the  N(2Z>)  4-  H,  reaction,  including 
experimental  study  of  NH(u,/)  nascent  rotational  states, 
which  could  not  be  probed  in  the  current  work.  Finally,  this 
investigation  points  up  the  need  for  an  accurate  determina¬ 
tion  of  the  H-N-H  potential  surface  in  order  to  predict  the 
reaction  dynamics. 

VIII.  SUMMARY 

Rate  parameters  for  NH(u)  creation  and  vibrational  re¬ 
laxation  by  Ar,  N2,  and  H:  have  been  measured.  The  reac- 
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tion  N(2D)  4-  H,-»NH(A><3)  -+-  H  forms  NH  in  excited 
vibrational  states.  Relaxation  rate  constants  k,.^ 

(M  =  N2,Ar,H2)  have  been  measured  for  the  first  time.  The 
relaxation  rate  constants  for  N2  and  Ar  are  small,  with  a  v 
dependence  intermediate  between  linear  and  quadratic;  up¬ 
per  limits  were  determined  for  rate  constants  k„=  (H:). 

The  effects  of  beam-created  species  on  NH(u)  quenching 
have  been  examined  through  both  detailed  reaction  model¬ 
ing  and  observation.  Experimental  conditions  were  chosen 
which  minimized  such  effects. 

The  nascent  NH(u)  vibrational  state  distribution  has 
been  measured  as  the  ratio  of  fitted  creation  rates  for  v  =  1- 
3,  also  for  the  first  time.  A  statistical  model  using  the  surpris- 
al  formalism  has  been  used  to  obtain  an  estimate  for  the 
relative  yield  of  NH(t'  =  0),  which  is  not  detected  in  our 
experiment.  The  nascent  product  distribution  is  more  highly 
excited  than  the  statistical  prediction,  with  much  less  prod¬ 
uct  in  f  =  0  than  predicted.  The  surprisal  analysis  is  consis¬ 
tent  with  the  reaction  N(2Z>)  +  H2  as  the  single  operative 
mechanism  for  NH(t’)  formation,  and  suggests  that  the  re¬ 
action  involves  a  direct  H  atom  abstraction  rather  than  a 
long-lived  intermediate. 
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